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I. INTRODUCTION 



Since the conjecture of Witten [1] that strange quark matter(SQM), would be more stable 
than the normal nuclear matter, much theoretical effort has been made on the investigation 
of its properties and apphcations [2-13]. Because of the well-known difficulty of QCD in the 
non-perturbative domain, many effective models reflecting the characteristics of the strong 
interaction are used to study the SQM. They include the MIT bag model [2-4], the quark- 
meson coupling(QMC) model [5], the Friedberg-Lee(FL) soliton bag model [6], the chiral 
SU(3) quark model [7], the quark mass density-dependent (QMDD) model [8-11] and the 
quark mass density- and temperature- dependent model [12, 13], etc. In this paper, we will 
focus our attention on the QMDD model. 

The QMDD model was first suggested by G. N. Fowler, S. Raha, and R. M. Weiner. 
According to this model, the masses oiu, d quarks and strange quarks (and the corresponding 
anti-quarks) are given by 

B - 
ruq^- — {i^u,d,u,d), (1) 

B 

ms,s = rriso + - — , (2) 

where niso is the current mass of the strange quark, B is the vacuum energy density inside the 
bag, and ub — ^{nu + nd + ris) is the baryon number density, with na, rig representing the 
density of u quark, d quark, and s quark. The basic hypothesis Eqs. (1) and (2) in QMDD 
model can easily be understand from the quark confinement mechanism. A confinement 
potential kr (r^) must be added to a quark system in the phenomenological effective models 
because the the perturbative QCD can not give us the confinement solutions of quarks. The 
confinement potential prevents the quark from going to infinite or to the very large regions. 
The large volume means that the density is small. This mechanism of confinement can be 
mimicked through the requirement that the mass of an isolated quark becomes infinitely 
large so that the vacuum is unable to support it [11, 12]. This is just the physical picture 
given by Eqs (1) and (2). In fact, this confinement mechanism is very similar to that of 
MIT bag model. But the advantage of QMDD model is that it does not need to introduce 
a quark confined boundary condition as that of the MIT bag model. 

Although the QMDD model is successful for describing the properties of SQM [8-11], 
but it is still an ideal quark gas model. Compared to the usual ideal quark gas model, the 
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basic improvement of the QMDD model is that the quark masses depend on density and 
the quark confinement mechanism is mimicked. Obviously, if we hope to investigate the 
physical properties of nucleons and hyperons by means of the QMDD model, the quark- 
quark interactions must be considered. Following the line of QMC model, we will introduce 
quark and scalar a field nonlinear coupling self-consistently to improve the QMDD model 
in this paper. As a first step, we ignore the s-quark and consider the coupling of u and d 
quarks to a non-linear scalar field only. 

We hope to emphasize there are two basic differences between our improved quark mass 
density- dependent (IQMDD) model and the usual QMC model suggested first by Guichon 
[14] and developed by Thomas, saito [15] and Jennings. Jin [16]. Firstly, instead of a MIT 
bag for nucleon in QMC model, we do not need a MIT bag for nucleon in IQMDD model. 
The constraint of MIT bag boundary condition disappears in our formulae because the quark 
confinement mechanism has been established in Eqs. (1) and (2). This is very important 
because it should provide a reasonable starting point for many-body calculations. Secondly, 
the interaction between quark and scalar meson is limited in the bag regions for QMC model, 
bur for IQMDD model, this interaction is extended to the whole free space. In fact, our 
model is similar to that of the F-L model, but instead of massless quarks in the F-L model, 
the masses of u and d quarks in our model are given by Eqs. (1) and (2). 

This paper is organized as follows. The main formulae of IQMDD will be given in the 
next section. The numerical results for the ground state and a number of physical quantities, 
namely, the root-mean-squared (rms) charge radius r^, the magnetic moment //p, the ratio 
between the axial- vector and the vector /3-decay coupling constants of the nucleon qa/ Qv 
will be presented in the third section. In the fourth section we will study the lowest one 
particle excited states with even or odd parity. The last section is a summary and discussion. 

II. THE IMPROVED QUARK MASS DENSITY- DEPENDENT MODEL 

We now briefiy outline the model and its pertinent features below. The Hamiltonian 
density of the IQMDD model reads: 

H = « • V + /3K + + ^TT^ + \{y(yf + U{a), (3) 
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where niq — B/Sub is the mass of u(d) quark, 'a and P are the standard Dirac matrices 
, represents the quark quantum field (color and fiavor indices suppressed) satisfying the 
canonical anticommutation relations: 

{^{r,t),^i;(r',t)}^5'(r-r'), (4) 

and / is the coupling constant between the quark field ip and the meson field a. The a field 
is considered independent of time and consequently the commutator: 

[7r{V),a{V')]=0, (5) 

where tt is conjugate field of the scalar meson field. Hence, a is treated as a classical field. 
In Eq. (3) U{a) is the self interaction potential for a field, which has a phenomenological 
form: 

^(-) = 1-^ + 1-^ + 1-^ + 5, (6) 

and 

cl > 3C2C4, (7) 

to ensure that the absolute minimum of U{a) is at cr = a^ac 7^ 0- The bag constant B is 
introduced in order that 

U{a^ac) = and U{0) = B. (8) 
Without any lose of generality, we may choose C3 < 0, and therefore a^ac > 0: 

(^vac = 7r-[-C3 + (C3' - |c2C4)'/']. (9) 

ZC4 6 

One can construct a Fock space of quark states and expand the operator ip in terms of 
annihilation and creation operators on this space with c-number spinor functions (pit^ , which 
satisfies the Dirac equation: 

[-7?- V + Pirn, + /a)]^W = ±e,9^W (10) 
z 

with superscripts ± denoting the positive and negative energy solutions, respectively. The 
spinor functions are normalized according to 
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^+^nd'r = 1. (11) 
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The total energy of the quark-scalar field system is given by 

E{a) = J2'- + I [^(V<7)^ + U{a)]d'r. (12) 



Minimum of E{a) occurs when a is the solution of 

-VV+^ = -/5^^„^„. (13) 

n 

The equations of the quark field for ground state and excited states will be given in Sec. 
3 and Sec. 4 respectively. 

III. THE GROUND STATE SOLUTION 

We discuss the ground state solution of the system now. Define (p as the lowest positive 
energy wave function. It can be expressed as: 

^ = { I Xm, (14) 



_ ,1 \ / . 

where a are the Pauli matrices, Xm — \ • The radial functions u{r) and v{r) 

\0 I I 1 

satisfy: 

^^-[e + m, + fairMr), (15) 

^ ^ -2^M + [6 - m, - fa{r)Hr), (16) 
ar r 

the normalized condition reads: 

POO 

Att [u'^{r) + v^{r)ydr = 1, (17) 
Jo 

and the equation of motion of a field Eq. (13) becomes: 

Equations (15), (16) and (18) can be solved with the boundary conditions: 

v{r = 0) = 0, u{r = oo) = 0, v{r = oo) = 0, a{r = oo) = ay, a'{r = 0) = (19) 



Before numerical calculation, we address the parameters of the model first. There are four 
free parameters, namely, C2, C3, C4, and / in IQMDD model. The parameters C2, C3 and C4 fix 
the interaction potential U, while / measures the coupling between the quark and the scalar 
field. The set of equations (15), (16) and (18) should be solved alternately until consistency 
was obtained using the iterative method [17]. Once the solutions of the above equations are 
obtained, one can calculate a number of physical quantities pertaining to the three-quark 
system which have been measured from experiments. Let r^, /ip, and Qa/ gv be, respectively, 
the rms charge radius, the magnetic moment, and the ratio between the axial-vector and 
the vector /3-decay coupling constants of the nucleon. They satisfy [18]: 

< >= j ip-^ipr^d^rj J cp-^cpd^ (20) 



^{Jfx ^^a^pdh)J j ip+ipd\ (21) 
and 

9a/ gv = ^J ^^a.ipd^r/ J ip+ipdh, (22) 
By using Eq. (14), we find: 

poo 

>=47r / {u^ + v'^ydr, (23) 
^0 



< rp^ >-- 



Afp = Y r^uvdr, (24) 



and 



gA/gv^'^j^ r^iu^-^v^r. (25) 

Since in the present form the model is fiavor independent, the charge radius of the neutron 
in the model is 

< rl >= 0, (26) 

and the neutron magnetic moment is 

2 

l^n^-^f^p, (27) 



6 



as given by the SU (6) algebra. Corrections to these relations arise only when QCD effects 
are included. 

For comparison let us remind ourselves that the experimental values of proton mass 
E = 4.69 fm~^, the proton magnetic moment n = 2.79 nuclear magnetons, Qa/Qv = 1-25, 
and the charge radius of proton Vp = 0.83 fm. To do numerical calculation, we follow Ref. 
[17, 18]. We fix Vp — 0.83 fm first. Once the value is chosen, there are only three free 
parameters. We have studied mostly two families of parameters for ground system state. 
These are the C2 — and the S^/^ = 145 MeV cases respectively. Each choice fixes one other 
parameter by a specific requirement for the shape of the potential. Therefore it is sufficient 
to label our results by the coupling constant / and constant C4. 

Our results for the two families (c2 = 0, B^^^ = 145 McV) are summarized in three tables. 
Through these tables, we use the units h — c — 1 and use fm as the fundamental unit for 
length. 

We consider the case B^^^ — 145 MeV first. The bag constant B corresponds to the 
difference between two minima of the U{a). In Table. 1 we fist the bag properties as a 
function of the parameter C4 for two values of the coupling constant /. The variation of bag 
properties with the coupling constant / for several values of C4 is given in Table. 2. Several 
features emerge from these calculations. Firstly, we note that an increase of the coupling 
constant produces a continuous change from a volume quark distribution for small / to a 
surface quark distribution for large /. This change is illustrated in Fig. 1 where we plot 
the quark density — for three values of /. It should be noted that shape of the sohton 
field does not change significantly with / for given C4. The variations of the quark charge 
density as a function of the radius with increasing C4 are plotted in Fig. 2 and Fig. 3. When 
C4 decreases, with a fixed value of /, the quark charge distribution — change slowly 
from surface to volume. This is illustrated in Fig. 2 and Fig. 3 for / = 75 and / = 200, 
respectively. The change from volume to surface quark charge density is also evident in the 
variations of the values of the magnetic moment //p and qa/ Qv- For instance from Table. 2, 
when C4 = 8 X 10^, the magnetic moment varies with increasing / from 1.87/xs to 2.61/xs, 
where /i^ is the Bohr magneton of the proton. Similarly, Qa/ Qv varies from 1.21 for small 
/ to 0.618 for large /. It is also found from the tables that for some C4 and / cases we can 
not found a solution for — 0.83 fm. 

The variation of the scalar field cr as a function of the radius for / = 50, C4 = 2 x 10^ 
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and B^/^ — 145 MeV is presented in Fig. 4. The value of a inside the hadron is very 
different from that of outside: inside a is less than zero, but outside a approaches (T^ac- The 
abrupt transition of scalar field through hadron surface will contribute to the total energy 
remarkably, as shown in Eqs. (11) and (12). 

TABLE 1. Variation of the properties with increasing parameters 



for the B^^^ = 145 MeV bag for the two values of /. 



/ 


C4 


e 


E 




9a/ gv 


30 


8 X 10^ 


1.54 


6.54 


2.39 


0.904 




1 X 10^ 


1.56 


6.56 


2.38 


0.929 




4 X 10^ 


1.74 


6.94 


2.25 


1.06 




8 X 10^ 


1.85 


7.26 


2.15 


1.11 


200 


8 X 10^ 




No solution 






1 X 10^ 




No solution 






4 X 10^ 


1.24 


5.23 


2.64 


0.603 




8 X 10^ 


1.22 


5.21 


2.63 


0.618 
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TABLE 2. Variation of the properties as a function of / for 
several values of C4 with B^^^ — 145 MeV. 





J 


C 


E 


H'P 


yAi yv 




1 

±0 




8 1 fi 
o.xu 


1 87 


1 91 

X .Z X 








7 9fi 


9 14 

Zi . X^ 


111 

X . X X 




■ifl 

uu 




fi 47 


9 3Q 


Q36 

U.Jl/OU 




1 


1 41 


5 Q7 


9 4Q 


7Q4 

U. / C/t: 




1 on 


1 94 


^ 70 

. 1 U 


9 ^'S 


790 

U . 1 zu 






1 99 


^ 91 


9 61 

Z/ . U X 


618 
u.uxo 


4 V 1 




9 1 1 

. X X 


8 01 
o.ux 


1 Q6 

X .C/U 


1 99 






1 74 

X . / ^ 


fi Q4 


9 95 

Z .ZO 


1 06 

X .uu 






1 4Q 


fi 98 


9 43 


868 
u.ouo 






1 3fi 


fi 07 

u.u / 


9 55 

Zi .00 


743 

U. / TlO 




1 Of) 


1 '^O 
X . ou 


^ 81 

.0 X 


9 57 

Zi .0 1 


680 
u .uou 






1 94 


^ 9'^ 

.ZiO 


9 6'? 

Z( .UO 


60"^ 
u .uuo 


1 V 10^ 


1 ^ 


1 Qfi 


7.62 


2.08 


1 1 7 

J — L / 




OKJ 


1 5fi 

X .OU 


6.56 


2.39 


Q9Q 






1 41 

X ."^tX 


6.09 


2.48 


u. / uo 




t 


x.oo 


5.77 


2.50 


n R77 




1 fin 




No solut 


ion 










No solution 




8 X 10^ 


15 


1.94 


7.61 


2.09 


1.16 




30 


1.54 


6.54 


2.39 


0.904 




50 




No solution 






75 




No solution 






100 




No solution 






200 




No solution 





Another family of parameters considered is characterized by C2 = 0. In this case U(o") 
has an inflection point at a = 0, and only one minimum. We vary C3, C4 and / subject to 
the bag size constraint Vp — 0.83 fm. Our numerical results for this family are summarized 
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in Table. 3. We find that when / increases, the character of the bag changes from volume 
confinement to surface confinement again. Similarly, we also find that for a given /, when 
C4 increases, the magnetic moment decreases. In summary, we find in both cases: 

(1) The quark energy e and the ratio Qa/qv decreases with increasing parameters / and 

C4. 

(2) On the other hand, the magnetic moment has just the opposite behavior: it is a 
decreasing function of the parameter C4 and an increasing function of the parameter /. 

(3) The total energy of the ground state decreases as a function of / and increases as a 
function of C4. The rise of the total energy E is caused by a contribution from 3e mainly. 

(4) An increase of the coupling constant produces a continuous change from a volume 
quark distribution to a surface quark distribution, and an decrease of C4 can reproduce 
same results. 
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TABLE 3. Variation of the properties as a function of / and C4 with C2—0. 





f 
J 


c 


E 


H-p 


ft A / rij T 

yAi yv 


1 y in^ 


1 ^ 


1 fi4 


7 nn 


9 31 


1 nn 

X .uu 




3n 


X .00 


fi 1 7 


9 "in 


n 7fiQ 

U. / Uc; 






1 9Q 


^ 77 


9 "17 
^ .u / 


n fifi7 

u.uu / 




/ 






9 fin 


n fi9i 




1 nn 


1 94 


5 43 


9 fil 

. U X 


n fini 

u.uux 








31 


9 fi3 


n 583 
yj .000 






1 fi4 


7 nn 


9 93 


1 n5 

X .\JO 






±.00 


fi 1 7 

U. X / 


9 48 


n 8in 

u.oxu 






1 9Q 


^ 77 


9 "17 
^ .u / 


n fiQ9 






1 9'i 


5 53 


9 fin 

.UU 


n fi35 




1 nn 


1 94 


5 43 


9 fi3 

Zi .UO 


n finp 

U.UUC/ 




1 "in 


1 99 


5 31 


9 fi3 


n 587 




1 ^ 


1 H'^ 
X . 00 


7 3fi 

1 .0\J 


9 Ifi 


1 1 9 

X . X ^ 




an 


1.48 


6.41 


2.46 


n.857 




sn 


1.34 


5.97 


2.55 


n.723 




75 


1.28 


5.71 


2.6n 


n.652 




inn 


1.25 


5.57 


2.63 


n.62n 




i5n 


1.24 


5.51 


2.61 


n.592 



IV. THE ONE PARTICLE EXCITED STATES 

The one particle excited states for F-L model had been calculated by Saly and Sundaresan 
[19]. In this section, we use their method to study the one particle excited state for IQMDD 
model. 

A. First excited state with even parity 

In this configuration, we shall consider two quarks to be in the ground state, while one 
quark will be placed in the first excited level ei. Then the system of equations to be solved 
is 
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^ = -[eo + m, + /<7(r)]^;(r), (28.a) 

^ = -2^ + [eo - m, - /<7(r))]i.(r), (28.6) 
diii(r) 



= -[el + mg + /a(r)]^;l(r), (28.c) 



^ = + [61 - m, - /a(r)]«i(r), (28.ci) 

ar r 

Att [u^{r) + v^{r)ydr = 1, (28.e) 



oo 



47r / [Mf(r) + i;f(r)]r^c^r = 1, (28./) 







^ + = + 2f[^?{r) r,\r)] + /[n?(r) - vl{r)]. {28.g) 

The set of equations (28.a)-(28.g) is to be solved with the boundary conditions: 



(29) 



v{r = 0) = 0, f i(r = 0) = 0, u{r = oo) = 0, v{r = oo) = 0, 
Vi{r — oo) = 0, Ui{r — oo) = 0, a{r — oo) = cry, a'{r = 0) = 0. 

B. Lowest energy with odd-parity state 

To obtain the lowest energy odd-parity state, we place two quarks in the ground state of 
even-parity and one quark in the lowest odd-parity state. In this case the system of equation 
is: 

^ = -[eo + m, + fa{rMr), (30.a) 
dv{r) _ _^v(r)_ ^ _ _ (30_^) 



dr r 
dvi{r) 



dr 



= [ei -ruq- fa{r)\ui{r), (30.c) 



du^{r) ^ _2«i(r) _ ^^^^ fa{r)\v,{rl (SO.d) 



dr r 

r"00 



POO 

An / [u^ir) + v^{r)ydr = 1, (30.e) 
Jo 

POO 

471 / ['Ui^(r) + vi^{r)ydr = 1, (30./) 
Jo 
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where Ui{r), Vi{r) are the components of wave function (/?: 



(30.<7) 



\ 





iviij-) 




The set of equations (30.a)-(30.g) is to be solved with the boundary conditions: 

v{r- = 0) = 0, Uiir = 0) = 0, u{r = oo) = 0, v{r = oo) = 0, 

Ui{r — oo) = 0, Vi{r — oo) = 0, a{r — oo) = cry, a'{r = 0) = 0. 

For simphcity, we consider S^/^ = 145 MeV, Vp — 0.83 fm case only. 

TABLE 4. Dependence of the solutions on C4 for /=30. 



(31) 



(32) 





Ground state 




Even-parity state 


odd-parity state 


C4 


e 


E 


eo 


ei 


E+ eo 


el E_ 


8 X 10^ 


1.85 


7.26 


2.02 


3.92 


9.11 1.87 


3.10 8.70 


6 X 10^ 


1.81 


7.16 


1.99 


3.99 


9.17 1.84 


3.11 8.71 


4 X 10^ 


1.76 


7.09 


1.95 


4.12 


9.29 1.80 


3.11 8.74 


2 X 10^ 


1.63 


6.72 


1.82 


4.16 


9.13 


No solution 


1 X 10^ 


1.56 


6.56 


1.72 


4.43 


9.36 


No solution 


8 X 10^ 


1.54 


6.54 


1.58 


4.41 


9.84 


No solution 




TABLE 5. 


De])eiideiiee of the solutions on / for r 


1=1 X 10^ 




Ground state 




Even-parity state 


odd-parity state 


/ 


e 


E 


eo 


ei 


E+ eo 


e[ E_ 


15 


2.09 


7.93 




No solut 


ion 


No solution 


20 


1.95 


7.53 


2.18 


3.80 


9.17 2.01 


3.21 8.82 


25 


1.85 


7.32 


2.05 


3.91 


9.21 1.88 


3.14 8.79 


30 


1.76 


7.09 


1.95 


4.12 


9.29 


No solution 


35 


1.66 


6.79 


1.80 


4.32 


9.29 


No solution 


40 


1.59 


6.59 


1.59 


4.28 


9.60 


No solution 


50 


1.49 


6.30 


1.48 


4.19 


9.63 


No solution 
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Our results are summarized in Tables. 4, 5 and Figs. 5, 6. In Table. 4,we fix = 0.83 
fm, 5^/^=145 MeV, / = 30 and show the dependence of ground state energy E, the first 
even-parity excited state energy E+, and the first odd-parity state energy E_ on C4. The 
dependence of E, E_|_, E_ on / for C4 = 4 x 10^ is shown in Table. 5. Remember that the 
experimental values of (uud) system are E = 4.69 fm~^, E_|_ = 7.35 fm~^ and E_ = 7.67 
fm~^ and compare our results with that given by Table. VI, VII of Ref. [19] for F-L model 
with massless quarks, we find that the agreement with experiments for IQMDD model is 
better than that of F-L model. 

Finally, we hope to point, as shown in Tables. 4 and 5, that the allowed range of 
parameters is large for the existence of the first even-parity excited state. The only restriction 
is that the height of the potential well fay must be greater than the quark eigenvalues Ci. 
But for odd-parity states, the allowed range of parameters restricts severely. We cannot find 
solution in many cases. To show this point more transparently, for odd-parity solutions, 
we fix Vp — 0.83 fm, B^^'^—lAb MeV, / = 30 and draw the curves of total quark density 
[2(m^ - v'^) -\- ui^ — vi'] vs. r and scalar field a vs. r for different C4 = 4 x 10^, 6 x 10^ and 
8 X 10^ in Fig. 5 respectively. We find no solutions exist for smaller values of C4. The same 
curves for fix C4 = 4 X 10^ but different / = 20 and 25 are shown in Fig. 6. We also find no 
solutions exist for higher values of /. 



V. SUMMARY AND DISCUSSION 



After introducing the quark and non-linear scalar field coupling, we suggest an improved 
quark mass density- dependent model. We obtain soliton solutions of ground state and 
excited states for the coupled equations for quark and scalar fields satisfying the required 
boundary conditions by numerical method. Since the most interesting question is whether 
the results of these calculations resemble the physics we are trying to describe, we concentrate 
on the dependence of the results on the phenomenological parameters introduced in this 
model. We present these dependence in the form of tables and figures for the ground state 
and low-lying excited states. The wave functions of quark are given. By using the wave 
function of ground state, we have calculated the rms charge radius, the magnetic moment 
and the ratio between the axial-vector and the vector /9-decay coupling constant of the 
nucleon and compared these values with experiment. We find the results given by IQMDD 
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model are in agreement with experiment. 

Since the boundary condition of MIT bag model has been given up in IQMDD model, 
the many-body calculation beyond mean field approximation can easily be carried out in 
this model. 

We note that the study of this paper is still limited at zero temperature and u, d quarks. 
Since the spontaneous breaking symmetry of U{a) will be restored at finite temperature, it 
is of interest to extend our discussion to finite temperature and study the effect of s quarks 
for hyperons. Working on this topic is in progress. 
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FIG. 2: Quark density — versus radius for C4 = 2 x 10^ and 8 x 10^, 5^/^=145 MeV and 
/=75. 
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FIG. 3: The same as Fig.2, but /=200. 
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FIG. 4: The soliton field a versus radius for /=50 and C4=2 x 10^ 5V4=i45 MeV 
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FIG. 5: Dependence of the odd-parity solutions on the parameter C4. Other parameters are kept 

at /=30, rp=0.83 fm for B^f^ = 145 MeV. No solutions exist for smaller values of C4. 
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FIG.6 



FIG. 6: Dependence of the odd-parity solutions on the parameter C4. Other parameters are kept 

at /=30, rp=0.83 fm for B^f^ = 145 MeV. No solutions exist for higher values of /. 
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